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Abstract 

The normalization of energy divergent Weber waves and finite energy Weber- 
Gauss beams is reported. The well-known Bessel and Mathieu waves are used 
to derive the integral relations between circular, elliptic, and parabolic waves 
and to present the Bessel and Mathieu wave decomposition of the Weber waves. 
The efficiency to approximate a Weber-Gauss beam as a finite superposition 
of Bessel-Gauss beams is also given. 
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1. Introduction 

Ideal propagation invariant scalar waves are the separable solution to the Helmholtz equa- 
tion with cylindrical symmetry. There exist four fundamental families related to the four 
cylindrical coordinate systems for which the reduced wave equation is separable: Plane 
waves for Cartesian symmetry, Bessel, Mathieu and Weber waves for circular-, elliptic- and 
parabolic-cylindrical symmetries, in that order [1, 2, 3, 4]. In optics, the interest in these 
families of structured scalar waves nowadays lies beyond their propagation invariance [5, 6]. 
It is common knowledge that they all carry a well defined linear momentum in the prop- 
agation direction. In addition, Plane waves carry a well defined linear momentum in the 
x and y directions, Bessel waves carry a well defined orbital angular momentum, Mathieu 
waves carry a well defined composition of orbital angular and x-linear momenta, while Weber 
waves carry a composition of orbital angular and y-\meai momenta [7, 8]. These dynamical 
properties are acquired by the corresponding beams and vector fields, both in the classical 
and quantum regime [9, 10, 11, 12], and can be used for manipulation of matter at the 
different scales [13]. 

These four propagation invariant families are exact solutions with divergent energies. By 
using the slowly varying envelope approximation (SVEA), these fundamental families have 
been shown to support finite energy solutions to Helmholtz equation [14]; e. g. Hermite- 
Gaussian beams are the most known solutions in the Cartesian coordinates. Finite energy 
beam families are also reported for the cylindrical symmetries, named as Bessel-, Mathieu- 
and Weber-Gauss beams, which are a closer description for the optical beams produced in the 
laboratory, e. g. those using holographic schemes [15] or optical resonators [16]. A theoretical 
description close to the experimental schemes of matter manipulation by these families of 
structured light requires the use of finite energy beams. In practice, the normalization of 
these structured light gives the requirement to determine the irradiance output needed from 
the laser sources to produce the desired interchange of mechanical variables. Furthermore, 
quantization of the corresponding finite energy fields requires the normalization of the scalar 
wave families. Currently, the general understanding of the parabolic waves, beams and fields 
includes their propagation characteristics and dynamical variables as well as holographic 
schemes for beam production [10, 12, 14, 15, 17]. But to one's surprise, the normalization 
for Weber-Gauss beams is missing in the literature. 

In this work, the link between Weber waves and the well-known Bessel and Mathieu 
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waves is provided by deriving the integral relations between parabolic, circular and elliptical 
waves following the phase space method proposed by Boyer, Kalnis and Miller [4, 8]. With 
the Bessel and Mathieu wave decompositions of Weber waves, we present the normalization 
of Weber beams under a series scheme [18]. The efficiency to approximate Weber-Gauss 
beams as a finite superposition of Bessel-Gauss beams is also given. The results found here 
should provide the necessary information for the quantization of the corresponding finite 
energy vector fields and for the understanding of any possible experimental demonstration 
of mechanical transfer involving such fields. 

2. Weber waves 

The wave equation in parabolic-cylindrical coordinates, u + t v = [2(x + % y)]^^ 2 \ accepts 
separable solutions of the form 
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These expressions, Weber waves, are Dirac delta normalized for a given frequency uo 

J d 3 ^Z^Zla = <W(7 - 7)5(5 - a). (3) 

The label set {p, k, 7, a} stands for the parity with respect to coordinate variables u and 
v, even or odd, wave number k = u/c, Euler angle corresponding to the decomposition 
of the wave vector in longitudinal and perpendicular components k = k± e± + k z e z = 
ksm r y e_i_ + A; cos 7 e z , and the real continuous eigenvalue, a, of the even, 

and odd functions, 

The definition of the hypergeometric 1F1 function follows Ref. [19]. All special functions 
notation and definitions will follow the latter reference. All calculations are based on the 
identities presented by Ref. [19, 20]. 
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Weber waves are eigenfunctions for the z-component of the linear momentum and the 
Poisson bracket of the z- and ^-component of the angular and linear momenta, in that order, 

Pz^ p,fc,7,a k z ^&p,k,"/,ai (8) 
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where the notation P = — «V and J = r x P has been used for linear and angular momentum. 
Thus, choosing a positive (negative) real eigenvalue a leads to horizontal parabolas opening 
to the left (right) with symmetry axis given by the x-axis. 

The Plane wave decomposition calculated for these waves is given by 
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where the angular spectra is written 
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Solutions without well defined parity can be constructed such that, 
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with angular spectra, 
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where the notation represents the Heaviside Theta function. The latter spectra are equiv- 
alent to those expressions presented in Ref. [4, 8]. 

Figure 1 shows a sampler of Weber waves for a given wave vector as well as positive and 
negative eigenvalue a. 
A. Bessel Decomposition 

Using the plane wave decomposition for Weber and Bessel waves, it is possible to show that 
the Bessel wave decomposition of a Weber wave is given by the superposition 
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Figure 1. Ideal Weber waves with wave number k = 2tt/X and z— component 7 = 0.31756 
eac 
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where the coefficients for the expansion are, 
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with auxiliary functions, 
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Notice that the complex conjugate of the auxiliary functions are given by C*(n, a) = 
C(n,—a) and f*(n,a) = f(n,—a), while the even and odd Bessel waves are defined in 
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the standard way, 

<<U 7l >> & *) = 0, z) ± *g,_ B (r, 0, *)], n = 0, 1, 2, 3 . . . (22) 

as functions of the Dirac delta normalized Bessel wave, 

* { 3, n (r, 0, z) = (2n)W>i n J n (k ± r)e< n * +k * z \ n = 0, ±1, ±2, . . . (23) 
with angular spectra given by the expression 

5. Mathieu Decomposition 

Eigenfunctions to the wave equation in elliptic-cylindrical coordinates, x+iy = f cosh[£+z7/], 
where / is half the interfocal distance of the coordinate system, are given by Mathieu waves 

^lUC V, z) = s e , n , q Ce m {^ q)ce m (v, q)e< k "~^\ (25) 
*$? 7 ,«& V, z) = s , n>q Se m (C, q)se m (v, q)e< k ^\ (26) 

where the parameter is defined as q = (fk±/2) 2 and the normalization coefficients for the 
four possible families, two parities (even or odd) and two periodicities (it or 2ir), are [21] 
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An interesting feature of Mathieu waves is that they take the form of ellipses (hyperbolas) 
as the relation a m (q) — 2q or b m (q) — 2q is positive (negative) with a m (q) the characteristic 
value for even Mathieu waves and b m (q) for odd. 

Mathieu waves have an angular spectra given by the expressions, 
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As the ordinary even or odd Mathieu functions are defined as a cosine or sine series, 

m 
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it is possible to relate the coefficients of the Mathieu decomposition with the aforementioned 
Bessel decomposition coefficients, 
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where the Mathieu decomposition of a Weber wave is given by 
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3. Weber beams 

A Weber beam is given by the SVEA solution to the Helmholtz equation [14], 

*W(r) = -j^e\ n «*)'e =^*Wa (M) . (38) 

The exponential part accounts for the Gaussian envelope with minimum waist wo and pa- 
rameter fj,(z) = 1 + iz/zr where the Rayleigh distance is given by zr = koj\j1 and the 
modified coordinates are defined as functions of x — x//i(z) and y = y//i(z). Please notice 
that the structure of the Weber waves is kept at the z = plane but it may change with 
propagation. 

A. Normalization. 

Normalization at the z = plane yields complex integrals wich can be avoided using a nor- 
malization scheme based on the Bessel decomposition [18]. Under this scheme, normalization 
for Weber beams yield 

oo 
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where the following notation has been used, oj = cj /2, 5 m , n stands for Kronecker delta, and 
I v for the z/th-order modified Bessel function of the first kind. 
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The series is convergent, in general it is possible to numerically argue that 

lim ^ < 1, (40) 

for any given value of a with 
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For the special case a = it is straightforward to show that 
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and the asymptotic limits [22] 
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Figure 2 shows the behaviour of the truncated normalization coefficient, ^j=o^' as a 
function of the total number of accounted terms, n, for even and odd Weber beams for a given 
wave vector and some random values of a. A thorough sampling in the range a G [—25, 25] 
showed that the truncated normalization coefficient was well stabilized at around the fiftieth 
term, at most. 
B. Generation. 

The holographic generation of a Weber wave or beam implies a great loss of irrandiance. An 
optical resonator with parabolic modes as output has not been reported yet to the knowledge 
of the author. Thus, a feasible scheme for obtaining a high irradiance Weber beam relies on 
the superposition of Bessel or Mathieu beams from optical resonators [16] . 
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Figure 2. (Color online) Behaviour of the truncated normalization coefficient Y^=o Sj f° r 
a given wave number k = 2tt/X and Euler angle 7 = 0.31756 rad leading to k z = 0.95 k 
for a) Even Weber fields ^^,7 a anc ^ ^) ^0^70=2 w ^h eigenvalues a = —11.2821 (circle), 
a = —2.2827 (square), a = (diamond), a = 10 (triangle), a = 23.7902 (upside-down 
triangle). 



In order to realize how many Bessel or Mathieu beams are required to reproduce a Weber 
beam, the behaviour of the coefficients for the Bessel or Mathieu decomposition has to be 
studied. Although the asymptotic limit for the modified Bessel function of the first kind, 
the Hypergeometric and the Gamma function are known [22, 23], it is quite complex to get 
an analytical asymptotic behaviour for the ratio between Bessel or Mathieu decomposition 
coefficients. 

For the Bessel decomposition, a thorough numerical survey for large values of n in the 
range a G [—25, 25] shows that 

lim 
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This can only be analytically asserted for a = 0. Figure 3 shows the numerical behaviour 
for some Bessel decomposition coefficients. It is possible to see that the main contributing 
terms are not the first terms of the series, with the exception of a = 0, as n increases the 
contribution of each term becomes important until it reaches a maximum, then it oscillates 
to stabilize and fulfil Eq.(49). 
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Figure 3. (Color online) Behaviour for the absolute value of the normalized Bessel decompo- 
sition coefficients of a) even and b) odd Weber waves with a given wave number k = 2tt/\, 
Euler angle 7 = 0.31756 rad leading to k z = 0.95 k, and eigenvalues a = —11.2821 (cir- 
cle), a = —2.2827 (square), a = (diamond), a = 10 (triangle), a = 23.7902 (upside-down 
triangle). 



4. Discussion 

The Dirac delta normalization of the even and odd ideal scalar parabolic waves given as 
hypergeometric functions, Weber waves, has been presented. The integral relations between 
these eigenfunctions of the Helmholtz equation with parabolic-cylindrical symmetry and 
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those with circular and elliptical-cylindrical symmetries, Bessel and Mathieu waves, have 
been shown and used to introduce the Bessel and Mathieu wave decomposition of Weber 
waves. 

A normalization for the finite energy Weber-Gauss beams was presented based on their 
Bessel-Gauss decomposition. It has been shown that it is not feasible to efficiently construct 
a Weber-Gauss beam through the finite superposition of just a few Bessel-Gauss beams. 

Finding a close analytical form for the normalization integral straight from the con- 
figuration or phase space representation of a Weber beam and the analysis pertaining the 
generation of Weber beams as the superposition of Mathieu beams is left as an open problem 
due to the complexity of the calculations involved. 

Acknowledgments 

The author is grateful to Prof. W. Miller for pointing to helpful references and thanks Prof. 
R. K. Lee, Prof. R. Jauregui and A. Stoffel for their useful comments. This work was 
supported by the National Tsing-Hua Univesity under contract No. 98N2309E1. 

References 

1. E. T. Whittaker and G. N. Watson, A Course of Modern Analysis (Cambridge Uni- 
versity Press, 1927). 

2. J. A. Stratton, Electromagnetic theory, International series in pure and applied physics 
(Read Books, 1941). 

3. P. M. Morse and H. Feshbach, Methods of theoretical physics, vol. 1, International series 
in pure and applied physics (McGraw-Hill, 1953). 

4. W. Miller, Symmetry and Separation of Variables, Encyclopedia of Mathematics and 
Its Applications (Cambridge University Press, 1984). 

5. J. Durnin, "Exact solutions for nondiffracting beams. I. the scalar theory," J. Opt. Soc. 
Am. A 4, 651 (1987). 

6. J. Durnin, J. J. Miceli, and J. H. Eberly, "Diffraction-free beams," Phys. Rev. Lett 58, 
1499 - 1501 (1987). 

7. A. A. Makarov, J. A. Smorodinsky, K. Valiev, and P. Winternitz, "A systematic search 
for nonrelativistic systems with dynamical symmetries," II Nuovo Cimento LII A, 4 
(1967). 



11 



8. C. P. Boyer, E. G. Kalnins, and W. M. Jr, "Symmetry and separation of variables for 
the Helmholtz and Laplace equations," Nagoya Math. J. 60, 35 - 80 (1976). 

9. R. Jauregui and S. Hacyan, "Quantum-mechanical properties of Bessel beams," Phys. 
Rev. A 71, 033411 (2005). 

10. K. Volke-Sepulveda and E. Ley-Koo, "General construction and connections of vector 
propagation invariant optical fields: TE and TM modes and polarization states," J. 
Opt. A: Pure Appl. Opt. 8, 867 - 877 (2006). 

11. B. M. Rodriguez-Lara and R. Jauregui, "Dynamical constants for electromagnetic fields 
with elliptic-cylindrical symmetry," Phys. Rev. A 78, 033813 (2008). 

12. B. M. Rodriguez-Lara and R. Jauregui, "Dynamical constants of structured photons 
with parabolic-cylindrical symmetry," Phys. Rev. A 79, 055806 (2009). 

13. D. L. Andrews, Structured Light and Its Applications: An Introduction to Phase- 
Structured Beams and Nanoscale Optical Forces (Elsevier, 2009). 

14. J. C. Gutierrez- Vega and M. A. Bandres, "Helmholtz-Gauss waves," J. Opt. Soc. Am. 
A 22, 289 (2005). 

15. C. Lopez-Mariscal, M. A. Bandres, and J. C. Gutierrez- Vega, "Observation of the 
experimental propagation properties of Helmholtz-Gauss beams," Optical Engineering 
456, 068001 (2006). 

16. M. B. Alvarez-Elizondo, R. Rodriguez-Masegosa, and J. C. Gutierrez- Vega, "Genera- 
tion of Mathieu-Gauss modes with an axicon-based laser resonator," Opt. Express 16, 
18770-18775 (2008). 

17. M. A. Bandres, J. C. Gutierrez- Vega, and S. Chavez-Cerda, "Parabolic nondiffracting 
optical wave fields," Opt. Lett. 29, 44 (2004). 

18. J. C. Gutierrez- Vega and M. A. Bandres, "Normalization of the Mathieu-Gauss optical 
beams," J. Opt. Soc. Am. A 24, 215 (2007). 

19. A. P. Prudnikov, J. A. Brychkov, and O. I. Marichev, Integrals and Series, vol. 3 
(Mockba, 1981). 

20. A. Erdelyi, ed., Higher Trascendental Functions, vol. 1 (McGraw-Hill, 1985). 

21. A. A. Inayat-Hussain, "Mathieu integral transforms," J. Math. Phys. 32, 669-675 
(1991). 

22. N. N. Lebedev, Special Functions and their Applications (Prentice-Hall, 1965). Trans- 
lator, R. A. Silverman. 



12 



M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Courier Dover 
Publications, 1970). 



13 



